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We show that the state with the highest known average two-particle 
von Neumann entanglement entropy proposed by Sudbery and one of the 
authors gives a local maximum of this entropy. We also show that this is 
not the case for an alternative highly entangled state proposed by Brown 
et al. 

1 Introduction 

The characterisation of multi-particle entanglement is a major open problem that 
is particularly significant for the study of quantum computation and many-body 
physics pp. One of the ways in which entanglement can be understood is by 
way of reference to a "maximally" entangled state. This target state can then 
be used, for example, to determine the largest rate at which it is possible to 
distill pure maximally entangled states from a supply of mixed states using only 
LOCC 0. An obvious condition for maximal entanglement in the case of pure 
states is that all one-qubit reduced density matrices are maximally mixed. For 
two- and three-qubit systems this leads to a unique state [3] (up to local unitary 
operations). However this is not true for a system with more than three qubits. 
For example, the states 
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have the property that all one-party reduced density matrices are maximally 
mixed; yet these two states are not locally equivalent. One may then ask which 
states also have maximally mixed two-party reduced density matrices. However, 
Sudbery and one of the authors (AH) have shown that it is not possible for all 
two-qubit reduced density matrices of a pure four-qubit state to be maximally 
mixed [4j. Nevertheless, they found a state which appears to maximise the aver- 
age von Neumann entropy of two-qubit reduced density matrices, which will be 
denoted by E2 in this paper. They showed that this state is a stationary point of 
the function E2, but it is not known whether this state indeed gives the maximum 
of E2. In this paper we show that this state gives at least a local maximum of E2. 

Given four qubits A, B, C and D the von Neumann entropy of the two-party 
reduced states is 

ExY = -tr{pxY log2 Pxy), (1) 

where pxY = iT^zw\4'){i^\, with W,X,Y,Z being a permutation of the four sys- 
tems. The average entropy of two-qubit reduced density matrices is defined by 

E2 = 77 (-^AB + Eac + Ead + Ebc + Ebd + Ecd) 
o 

= '^^^AB + Eac + Ead)- 

This quantity can naturally be taken as a measure of the entanglement between 
two pairs in pure four-qubit states. 

There are other approaches to quantifying multi-partite entanglement that are 
applicable to a general mixed state. A measure for a general composite system has 
been introduced by Yukalov [5J using the ratio of norms of an entangling operator 
and of a disentangling operator in the relevant disentangled Hilbert space. The 
entanglement of a four-qubit system can also be studied using the entropy of the 
reduced three-particle system and the strong subaddativity inequality [6]. Brown 
et al. have considered the partial transpose with respect to all possible partitions 
of the state [7]. These measures are more general than E2 in so much as they 
are applicable to mixed states. In this paper, we only consider pure states and, 
therefore, the average entropy is a suitable measure of entanglement. 

The four-qubit state proposed in 

IM4) = ^[|0011) + |1100) + cj(|0101) + |1010)) +U(|0110) + |1001))], 

where uj = exp(2z7r/3) is a third root of unity, has the highest known average two- 
qubit bipartite entanglement [HEj. It is also an element of the orbit of SLOCC 
operations that has maximal four-partite entanglement [H]. Together with its 
complex conjugate, it also provides a basis for the space of singlets contained in 
a four-qubit Hilbert space [9j. 

The entropy Exy measures the entanglement between systems XY and WZ, 
i.e. it measures entanglement between pairs. The entropy of the one party 
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reduced density matrices measure the entanglement between individual systems 
and the rest of the state. Another way in which entanglement could manifest itself 
in a four qubit system is the entanglement between any two individual systems. 
For example, the entanglement between X and Y is measured by regarding pxY 
as a (mixed) state in its own right. There are various bipartite entanglement 
measures for mixed states. One that is commonly used for multipartite states 
is the concurrence, since for this measure the CKW inequality gives an upper 
bound on the bipartite entanglement in terms of one party entanglement pUlfTT] . 
Namely, 

where Ca{bcd) denotes the concurrence across the partition A : BCD i.e. con- 
sidering the qubits BCD as one 8 dimensional qudit. Since pabcd is pure, 
Ca{bcd) = 2^/det{pA). 

It is interesting, to note that IM4) contains no bipartite entanglement, that 
is. Pay is separable! for all Y. Moreover, the states pAv are on the boundary of 
separable states. This can be seen by writing pay in the form 

fl - -) 1 

PAY = ^^I+-|$-)($-|, 

where |$_) = ^(|01) - |10)), i.e. as a Werner state with x = ^ [HI US]. There- 
fore, the l.h.s. of eq. ([2]) is zero while the r.h.s. is one. This can be thought of as 
saying that A shares none of it's entanglement with it's neighbours i.e. the en- 
tanglement between AY and WZ that cannot be regarded as the entanglements 
of A with W or Z separately is maximised by IM4). 

The authors of [Ij conjectured that IM4) gives a maximum of the average 
two-qubit bipartite entanglement. By considering the first variation, 6Exy, they 
were able to show that E2 is stationary at IM4). In this paper, we will consider 
the second variation of the average entropy and demonstrate that the state IM4) 
indeed gives a local maximum. We will also consider another highly entangled 
state proposed in |7j and show that this state is in fact not a stationary point 
of our measure, thus illustrating that maximal entanglement is dependent on the 
measure used. 



2 Varying the Entropy 

We consider variations of the state near IM4) with the varied states characterised 
by several small parameters. In general, suppose we vary a four-qubit state lip) 
in second-order approximation as 

^ Pay is a mixed state, and so separability here means that it can be writen as the convex 
sum of uncntanglcd pure states. 
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where \5ip) and are of first order and of second order, respectively, in the 

small parameters characterising the variations. From the normalisation condition 

{ip'\ip') = {iplip) = 1, we obtain 

2Re{6x/j\^) = 0, (3) 

at first order and 

2Re{5'^^p\^) + {6^p\6^|J) = 0, (4) 

at second order. 

In order to find the first and second variations of the entropy, we need the 
following lemma. 

Lemma 2.1. For any function f{A) of a matrix A that can be written as a 
power series, consider variations of A with small parameters, A ^-^ A + 6A + 6'^A, 
to second order in these parameters. Then to second order the corresponding 
variation tr[/(y4)] t-^ tr[/(y4)] + 5ti[f{A)] + 5hi[f{A)] is given by 

6tT[fiA)]=tT[6A-f'iA)], 
6hT[f{A)] = tT[6'A . f'iA) + ^6A ■ Sf{A)]. 

Proof. Since f{A) can be written as a power series, it is enough to show these 
formulas for f{A) = A". We obtain first-order terms in the variation of A" by 
replacing one of the A's by 6 A. Thus, StiA^ = ntT[6A ■ A"'^^]. This proves the 
first formula. We get second-order terms by replacing one of A's by a 6"^ A, or by 
replacing two A's by two 6A^s. Thus, 

n-2 
fc=0 

= tr [S^A ■ nA"-^] + ^tr [6A ■ . 
(One can perhaps convince oneself of the need for the factor 1/2 in the second term 



by noting that the number of ways to replace two A's by two 6A^s is n{n — l)/2.) 
This proves the second formula. □ 

Letting A = 1 — pxy and f{A) = {1 — A) log(l — A) = Pxy^o&Pxy in this 
Lemma, and using the normalisation conditions tr[^pxy] = tr [5^pxy] = 0, we 
obtain the first- and second-order variations of Exv defined by as 

SExY = -j^^tr[5pxy logpxy], (5) 
S'^ExY = -r^tr[5Vxy log Pxy + ^^Pxy 5 log Pxy], (6) 
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where 



PXY = tYzw{\i^){'ip\), 
6pxY = tTzw{\S^){^\ + \^){Si^\), 

3 An alternative highly entangled state 

In [7j, Brown et al. have considered multi-partite entanglement across all possible 
partitions of a state and calculated the sum of all negative eigenvalues when the 
partial-transpose function is applied and sought to maximise this over all possible 
states. In the four-qubit case, their numerical search found the state 

1^4) = ^(lOOOO) + I + oil) + 11101) + I - 110)), 

where |+) = -^(|0) + |1)) and |— ) = :^(|0) — |1)). Here we consider whether the 
state \1lj4) gives a local maximum of E2 as a measure of entanglement. 

We start by noting that p^c is maximally mixed and that the eigenvalues of 
Pab and pad are both equal to 

2 + V2 2 + V2 2-^2 2-V2I 
8 ' 8 ' 8 ' 8 J ■ 

Comparing the average two-qubit entanglement of this state with that of IM4), 
we see that the state 1-04) has lower two-party average von Neumann entangle- 
ment [7]. Namely, 

E2{\^P,)) = I - ^ log, (3 + v^) ^ 1.7426, 
E2(|M4)) = l + ^log22^ 1.7925. 

We will now show that the state in fact does not give a local maximum of 
E2 by demonstrating that it is not a stationary point. The first-order variation 
of E2 is given by 

35^2 = 6Eab + SEac + 5Ead. 

Since pac is maximally mixed, we have SEac = by ([5]) and the first-order 
normalisation condition. In order to diagonalise pab and pad "we use the bases 
{|^±)) k±)} and {\w±), \x±)}, respectively, where 

\u±) = (v^T 1)|10) ± |00), \v^) = (v^± 1)|11) T |01), 

\w±) = (v^T 1)|10) T |00), |x±) = (v^± 1)|11) ± |01). 
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Therefore we can now write 



pAB — Pad 



2 + ^2 



2 + V2 



2- V2. 



Hence, by ([5]) and the normahsation condition, 

5E2 oc ii[5pAB{\u+){u+\ + \v+){v+\) + 6pad{\w+){w+\ + 

= {u+\5pab\u+) + (f+|<5pAB|f+) + (w+|5p^d|w+) + {x+\5pad\x+). 

Consider a variation of the form 

1^^) = a|0011) +/?|1011). 

Then the normahsation condition requires that Re{a + /5} = 0. We find 

5pAB = 4^(|00) + |10))(a(00| + p{10\) + 4^(a|00) + /5|10))((00| + (10|), 



and hence 



^/2 V2 

{u+\Spab\u+) =2Re \a + P{V2-l] 
{v+\Spab\v+) = 0. 



Similarly, 



{w+\6pad\'w+) = 0, 

{x+\6pad\x+) = 2Re \a{V2 + 1) + (3{2V2 + 3) 
Therefore, for the variations considered here we have 



a{2 + V2) + /3{2 + 3V2 



6E2 oc Re 



Hence, by putting a = —(3 = e for a small £ G M — note that the normalisation 
condition Re(a + /?) = is satisfied — we have 6E2 7^ 0. Therefore the state \ip4) 
cannot give a local maximum of E2. 
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4 The Second-order variations 



We now return to the state IM4) and show that it gives a local maximum of E2, 
that is, S^Eab + S^Eac + S'^Ead < 0. Let us write S'^pxY = i^xv + o'xy, where 

axY = tTwz{\S^p){Si^\), 

Then, 

-log2'^ S'^Eay = ^tr(KAy logpAy + cTAy log p^y + ^^p^y^ log p^y), (7) 
y y 

where Y = B,C and D. Our task now is to show that the right-hand side of this 
equation is positive definite for all nontrivial variations of the state {ip) = IM4) 
satisfying the normalisation conditions and (jl]). We will deal with each term 
in ([7]) separately. 

4.1 The first two terms in the expansion 

We note that logp^y = log3- |$-)($-| -log6-I, where |$_) = ^TjdlO) - |01)), 
for all Y, and that 

tr(/€AB|$_)($_|) = Re{{6^^\^) - {6^^m, 

tr(/€Az?|$-)($-|) = Re((5'^|^) 

where \ip) is the complex conjugate of lip) = IM4) in the computation basis, as was 
shown in [4] in the context of first-order variation. Hence, using the normalisation 
condition (jl]), we find the first term of ([7]) as 

^tr(KAy log pAy) = 31og2\/3- {SiplSilj). 

Y 

We now consider the second term in ([71). We have 

tr (cTxy log pxy) = tr(crxy log3 • |<l>_)(<l>_| - crxy log 6 • I) 

= logs ■ tT[itTwzm{6M^-){^-\] - log6 ■ {5^\6iIj) 
1 

= logs 5^ |(($_Uy(zjVz)|5^)|' - log6 ■ {6tP\6ij). 

i,j=0 

Hence, 

^tr((TAy log pAy) = -Slog6 • {dipldip) + Fab + Fac + Fad, 

Y 
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where 

1 

FxY = 

i,j=0 

Thus 

J2 [(kay + (Jay) log pay] = -Slog Vs ■ {Si/j\5ij) + Fab + Fac + Fad 

Y 

because Fab + Fac + Fad > 0. This motivates us to define 

P = J2iT[SpAYSlogpAY] - 31og3 • {5ij\5ij) . (8) 

Y 

Then, if P > for aU nontrivial variations, 5^i?2 is negative definite and the 
state {M/i) gives a local maximum of We will show this fact with a certain 
convenient parametrisation of variations. 

4.2 The third term in the expansion 

The following lemma will be useful in analysing the variation t^logpxy- 

Lemma 4.1. Provided that the eigenvalues of A are positive and less than 1, we 
have ^ 

5\ogA^ f [I-t{I- A)]-^5A[I-t{I- A)]-^dt. 
Jo 

Proof. We expand log74 as 

log A = log(I -{l-A)) = -J2 -■ 



n=l 



Then 



(I - .4)"' .5.4(1 - .4)"' 



OO oo 



n' + m + 1 

n'=0 m=0 



where we have let n' = n — m — 1. Noting the elementary integral, 



^0 



n' + m + 1 ' 
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we find 

oo oo 



Y.^t{^ " A)r6A[t(I - A)fdt 

- n'=0 m=0 

oo oo 

m=0 n'=0 

as required. □ 

We use this lemma with A = p, where p is a density matrix. Since p is 
Hermitian, we can choose a basis in which p = diag{Ai, . . . , A„}. We can apply 
this lemma if < Aj < 1 for all i. Note that the state IM4) has this property. If 
Aj 7^ Xj, we have 



^^^°^^)-=/ i-ta-A/ ^' i-ta-A,) ^^ 

^ ^ X Y 







1 - t(l - Ai) 1 - t(l - Xj 



6pijdt, 



where X = j-^y Y = j-^- Hence, for Aj 7^ Aj, 



If Aj = Aj, then 



(.Slogp),, = j-^^ log (^y^ ) Spij. (9) 



A,; 



(10) 



(This formula can also be obtained by letting Xj — > Aj in ([9]).) We now apply 
these formulas to the variation 6 log pxY- 

To ease the notation we let 6pAB = ('^ij^)) ^Pac = ('^Ij'*) s-iid ^Pad = (c^l^'')- 
We will use the basis S = {|0), |/), |+), |-)} = {|00), |11), |$+), |<l>_)}, where 
1$^) = 2-i/2(|io) - |01)). Thus, for example, agj = {0\Spab\I)- Since Ai = 
A2 = A3 = 1/6 and A4 = 1/2 for each pay, we have, by applying and (fTUl) . 

ti{5pAB5 log Pab) = "Y^i: 



HjPij I ; 



tr((5p^c"5 log Pac) = 

ij 

ti{5 pad5 log Pad) = y^^K, 



u \^tj I ' 



«j I I ' 
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where 



/ 6 



6 
6 
6 

3 log 3 



6 
6 
6 

3 log 3 



31og3 \ 
3 log 3 
3 log 3 

2 / 



K 



6 
6 

V 31og3 



In order to proceed further, we need to explicitly parametrise the variations 
1 5^). Thus, we write 



\5^) = eoooo|0000) + enii|llll) 

+eooii|0011) + enoolllOO) + ^(eioio|1010) + eoioi|0101)) 
+^2(eiooi|1001) + eoiio|0110)) 

+eoiii|0111) + eioiillOll) + eiioi|1101) + emo|1110) 
+2(11000) + |0100) + |0010) + |0001)) . 



Note that we have included only one term with "three O's" out of four possible 
terms. This is because all other terms can be eliminated by local unitary trans- 
formations to first order. We derive additional constraints on the variations by 
noting that if the effect of the variation is to change the relative phase in any one 
of the qubits, then our new state lifj') is locally equivalent to IM4). Let us write 



The first-order variation in the relative phase within the first qubit results in the 
change in yi + y2 + Vs- Similarly, the phase variations in the second and third 
qubits change the values of —yi +1/2 + 2/3 and 2/1 — 2/2 + 2/35 respectively. Hence, 
for any variation, we can always find an equivalent variation satisfying 



by adjusting these phases. In the same way we find that an overall change of 
phase, IM4) e*^|M4), can be used to impose the condition 



Thus, we have 21 real parameters (after imposing the normalisation condition) 
in our space of variations. Since the dimensionality of the space of locally in- 
equivalent states is 18 (see, e.g. (TU]) in a neighbourhood of a generic state, three 
dimensions are redundant. This discrepancy is due to the fact that the state 
IM4) remains unchanged if all qubits are transformed by the same SU{2) ma- 
trix: thus, the dimensionality of the orbit of the local unitary transformations of 
the state IM4) is 10, which is smaller than the dimensionality of this orbit for a 



Cooii — ^1100 = Xi + iyi, eiioo + eooii = Xi + iYi, 
coioi ~ ^1010 = X2 + iy2, cioio + ^0101 = X2 + iY2, 

COIIO ~ ^1001 = X3 + Z2/3, eiooi + Cono = -^3 + ^^3- 



2/1 = 2/2 = 2/3 = 



(11) 



Y, + Y2 + Y, = 0. 



(12) 



10 



generic state by 3. We have eliminated all variations that reduce to infinitesimal 
local unitary transformations, but the set of physically equivalent variations is 
generically three dimensional. It would be possible to eliminate this redundancy 
by using canonical forms [Tl[ [15] though we have chosen not to do so. 

With the condition (ITTll imposed, the quantity P defined by ([8]) can be written 

as 

P = P^ + P^ + P^ + P^^ 

where 

3 

Pi = 12j]|aSr-31og3(|eoooor+|eiiiir), (13) 



3 



P, = 5^[l2(|aKr + |aSr) + 61og3(|afJr + 



Q = l 



-31og3(|eoiiir + leioiil + |eiioi| + |eiiio| + M ) , (14) 

3 o 

Pa = 61og3^|aJ'lp--log3(x? + x^ + x2), (15) 

a=l 

A = E[6(|aSP+l4fr + |anr) + 2|aL"2r]-^log3X:(X^ + 

a=l a=l 

(16) 

We will show that (i) Pi > if either eoooo or eim is nonzero, (ii) P2 > if any 
of eoiii, eioii, eiioi, emo or z is nonzero, (iii) P3 > if any of x^'s is nonzero and 
(iv) P4 > if any or X^'s or Y^s is nonzero. This will imply that P is positive 
definite. 

4.2.1 Posit ivity of Pi 

The only terms relevant here 's. These are given by 

VeagJ = {0\6pab\I) = mil) + (00|#) = eTTTT + eoooo, 
\fQa^o] = {0\6pAc\I) = uj{6iIj\II) + uJ{00\6ip) = t^elUT + Ueoooo, 
A/6a^| = {0\6pad\I) = Lj{6tlj\II) + uj{00\6iIj) = cjemi + t^eoooo- 



Thus we have 



3 ^ 

ko/p = ^(koooop + kiiiin 



a=l 



Hence by ([T3]) 

Pi = (6 - 31og3)(|eoooo|^ + |eiiii| 

which is positive if either eoooo or eim is nonzero. 
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4.2.2 Posit ivity of P2 

We find the relevant a-^'s as 



V^afl = --^(eiioi + eiiio) + ^2 2;, 



1+ 
(1) 



i (eiiio — ciioi); 



Hence, 





3, 






f + 




! 3 eiioi 


— ClllO 


f + 




; 3 eioii 


~ Coin 



12|(ij^| -|- 6 log 3 1 cij'_ I — i^iiui I "-iiiu I f-, ^"^6 '-'I'-iiui ^iiiui ) 

12|ag|p + 61og3|agl|2 = \e,on + ' 

The corresponding quantities involving a[f^s and a[f can be obtained similary 
as follows: 

12|ag|2 + 61og3|afi|2 
12|ag||2 + 61og3|ag^|2 

12|4^_J|2 + 61og3|af2|^ 
12|ag|p + 61og3|ag^|2 



,2 3 
1+2 


log 


;3 eioii 


— ClllO 


f + i 

' 2 


log 


;3 ciioi 


~ ^0111 



,2 3 
■1+2 


log 


; 3 eioii 


— eiioi 


' 2 


log 


; 3 eiiio 


~ ^0111 



From these equations we find 
1 
2 



P2 = ( 1 - ■^log3 ) (|eoiii + eioiil^ + |eoiii + eiioil^ + |eoiii + eiiioi 



+ ( 1 - i log3 ) (|eioii + eiioil^ + |eioii + emol'' + |eiioi + emol 



+ log3(|eioii — Ciioip + |eioii — emop + |eiioi — ei 



110 



+24 (^l-ilog3j \z\'. 

It is clear that the right-hand side is positive unless eom, eion, enoi, emo and z 
all vanish H 

^The expression for P2 is not symmetric under permutations of four qubits involving the first 
qubit. However, (5^ £'2 itself is symmetric under such permutations thanks to the contribution 
from Fab + Fac + Fad (which we have discarded because it is positive definite) . This must 
be the case because the average two-partite von Neumann entanglement entropy E2 has this 
symmetry. 
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4.2.3 Positivity of P3 

We have 

2\/6a^l = -2(x2 + X3) - V3i{x2 - X3), 

(2) a) 

and a\__ and a\__ are obtained from this by cychc permutations 2 — 3 ^ 1 and 
3 ^ 2 — *^ 1, respectively. Hence, 

Thus, P3 given by (|T^ is 
lo 3 

-P3 = —^{^^1 + ^^2 + 43^3 + X2X3 + X3X1 + X1X2) , 



which is positive if Xi, X2 or 2:3 is nonzero. 

4.2.4 Positivity of P4 

We have 



Hence 
Similarly, 



^^^00 ~ ^0011 + ^0011 — Xi + xi, 
VQa\j = eiioo + eiioo = -^1 — xi. 

6|agiP + 6|4yp = 2(X2 + a;?). 

6|agi|2 + 6|af;p = 2(X2 + x2), 
6|agir + 6|af;p = 2(X| + x^). 



Thus, 



a=l a=l 

The remaining 'diagonal terms' are 

-VQa^^l_ = Re{u;(eioio + eoioi) + ^^(eoiio + eiooi)}, 
-V2a^}l = Im{t^(eioio + eoioi) - ^^'(eoiio + eiooi)}- 

The coefficients a^^\_ and a*^! {0,++ and a^l) are obtained from the expressions 
for and a^}}_ by interchanging the second and third (fourth) qubits. Remem- 
bering the definitions eioio + eoioi = X2 + iY2 and eiooi + eono = X3 + iY^, we 
obtain 

6|aJ|p + 2\a^l\' = X| + + Xj + + 2(X2X3 - Y2Y,). 
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We find similarly 

6|aJiP + 2|aL'lp = Xl + + Xl + + 2{X,X, -Y,Y,), 
Q\at\\^ + 2|aL'l I' = Xl + Y^ + Xl + if + 2(XiX2 - Y1Y2). 

By combining these formulas and (JTTl) with the definition (fT6|) of P4 we have 
P4 > (^4-^log3^ (X2 + x2 + X|) + 2(X2X3 + X3Xi + XiX2) 
+ (^2 - ^ log 3^ {Y^ + + Fg^) - 2(F2>^3 + Y^Yr + Y^Y^). 
Finally, the use of the condition Fi + F2 + ^3 = leads to 

P4 > [(X2+X3)2 + (X3 + Xi)2 + (Xi+X2)2] 

+ (^2 - ^ log 3^ {Xl + X2 + Xl) + 3 (2 - log 3) (Y^^ + y^Y^ + F^')- 

Thus, P4 > unless Xq,'s and y^'s all vanish. This completes the proof that (5^i?2 
at IM4) is negative definiteEl Hence, the state IM4) indeed gives a local maximum 
of the average two-partite von Neumann entanglement entropy. 
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